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ABSTRACT:  Sarin  and  soman  are  warfare  nerve  agents  that  represent  some  of 
the  most  toxic  compounds  ever  synthesized.  The  extreme  risk  in  handling  such 
molecules  has,  until  now,  precluded  detailed  research  into  the  surface  chemistry  of 
agents.  We  have  developed  a  surface  science  approach  to  explore  the  fundamental 
nature  of  hydrogen  bonding  forces  between  these  agents  and  a  hydroxylated 
surface.  Infrared  spectroscopy  revealed  that  both  agents  adsorb  to  amorphous 
silica  through  the  formation  of  surprisingly  strong  hydrogen-bonding  interactions 
with  primarily  isolated  sOanol  groups  (SiOH).  Comparisons  with  previous  -3 
theoretical  results  reveal  that  this  bonding  occurs  almost  exclusively  through  the 
phosphoryl  oxygen  (P=0)  of  the  agent.  Temperature-programmed  desorption 
experiments  determined  that  the  activation  energy  for  hydrogen  bond  rupture  and 
desorption  of  sarin  and  soman  was  50  ±  2  and  52  ±  2  kj/mol,  respectively. 
Together  with  results  from  previous  studies  involving  other  phosphoryl-containing 
molecules,  we  have  constructed  a  detailed  understanding  of  the  structure- 
function  relationship  for  nerve  agent  hydrogen  bonding  at  the  gas— surface  interface. 

SECTION:  Surfaces,  Interfaces,  Porous  Materials,  and  Catalysis 
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The  threat  of  chemical  warfare  agents  (CWAs),  particularly 
the  nerve  agents  sarin  (RS-propan-2-yl  methylphosphono- 
fluoridate,  GB)  and  soman  (3,3-dimethylbutan-2-yl  methyl- 
phosphonofluoridate,  GD),  to  the  global  community  has 
motivated  researchers  to  explore  new  methods  for  the  detection 
of  agents,  the  destruction  of  existing  stockpiles,  and  the 
decontamination  of  affected  areas,  with  the  ultimate  goal  of 
providing  effective  protection  of  military  and  civihan  personnel. 
Critical  to  the  development  of  strategies  for  combating  the 
threat  of  CWAs  is  a  detailed  understanding  of  agent— surface 
chemistry.  However,  scientists  have  largely  been  restricted  to 
conducting  experiments  with  less  toxic  analogues  (often 
referred  to  as  simulants)  of  the  actual  agents.  These  simulants 
(Eigure  l)  contain  the  central  phosphoryl  ester  of  the  actual 
agents,  but  they  lack  the  structure  and  key  chemical 
functionality  of  the  agents,  which  makes  the  prediction  of 
agent— surface  chemistry  tenuous.  Eor  example,  recent  work  has 
shown  that  dimethyl  methylphosphonate  (DMMP),  the  most 
commonly  employed  simulant,  binds  to  the  surface  of 
hydrophilic  silica  much  more  strongly  than  dimethyl  chlor- 
ophosphate  (DMCP),^  which  implies  that  a  phosphorus- 
halogen  bond,  such  as  that  found  in  sarin,  may  affect  the  surface 
chemistry  of  the  agent  in  ways  that  cannot  be  predicted  by 
studying  the  simulant  alone.  In  fact,  calculations  suggest  that 
the  electron-withdrawing  character  of  the  fluorine  group  in 
sarin  has  a  major  impact  on  its  interfacial  hydrogen-bonding 
strength.^  We  have  therefore  developed  an  ultrahigh  vacuum 
(UHV)  surface  science  approach  to  the  study  of  agent— surface 
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Figure  1.  The  structures  of  the  CWAs  sarin  (red)  and  soman  (blue) 
and  five  of  the  most  common  simulants  (black)  used  to  help  predict 
agent  chemistry. 


chemistry.  Here^  we  present  direct  studies  of  the  mechanisms 
and  energetics  that  govern  hydrogen-bonding  forces  when  sarin 
and  soman  adsorb  to  one  of  the  most  abundant  materials  found 
on  Earth;  silica.  In  addition  to  the  relevance  of  sihca  surface 
chemistry  to  many  areas  of  environmental  science^  the  material 
serves  as  an  excellent  test  system  for  exploration  of  the 
fundamental  nature  of  molecule— surface  hydrogen  bond  forces^ 
which  play  a  critical  role  throughout  chemistry  and  control  the 
uptake  and  residence  times  of  molecules  on  the  surface  of  a 
variety  of  materials  from  coatings  to  oxides. 
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Experiments  were  performed  in  an  UHV  surface  science 
instrument  recently  developed  for  the  study  of  CWAs  and 
other  highly  toxic  compounds  (see  the  Experimental  Methods 
section).  Located  in  a  laboratory  at  the  Aberdeen  Proving 
Grounds,  the  instrument  (Figure  5)  enables  measurements  of 
surface  adsorption  through  infrared  (IR)  spectroscopy  and  X- 
ray  photoelectron  spectroscopy  (XPS),  while  gas-phase 
products  are  monitored  via  a  quadrupole  mass  spectrometer.^ 
Particulate  silica  surface  samples  were  prepared  by  dispersing 
silica  (200  m^/g,  Aerosil  fumed  silica  with  a  12  nm  average 
particle  diameter)  from  a  solution  of  isopropanol  onto  a 
reflective  Au  substrate,  which  served  as  an  inert  support  and 
reflective  mirror  for  performing  IR  spectroscopy  of  the  silica 
during  agent  exposure.  Once  prepared,  the  surface  samples 
were  installed  via  a  load-lock  mechanism  that  couples  the  main 
chamber  sample  holder  to  a  small  transfer  chamber  located 
within  the  confines  of  a  CWA-certified  surety  fume  hood. 
Within  the  main  chamber,  the  sample  was  mounted  on  a 
molybdenum  platen,  which  can  be  heated  via  a  tungsten 
filament  and  cooled  with  hquid  nitrogen.  Prior  to  each 
experiment,  the  sample  was  cleaned  by  thermal  annealing  to 
700  K  before  cooling  the  sample  to  200  K,  where  dosing  was 
conducted.  IR  and  XPS  measurements  were  used  to  verify  that 
the  surface  of  the  sihca  was  clean  to  less  than  0.1%  coverage  of 
organic  contaminants.  In  addition  to  sample  cleaning,  heating 
to  700  K  produces  a  partially  dehydroxylated  surface.  The  OH 
coverage  for  our  sample  is  estimated  to  be  2  OH/ nm^,  which  is 
typical  for  particulate  silica  under  these  sample  preparation 
conditions.^  ^  Chemical  Agent  Standard  Analytical  Reference 
Material  (CASARM)-grade  sarin  (98.7%  by  NMR)  and  soman 
(96.0%  by  NMR)  were  used  for  aU  experiments.  (CAUTION: 
Experiments  performed  with  ultratoxic  CWAs  require  extreme 
care,  highly  trained  operators,  and  a  secure  government- 
regulated  facility.)  Prior  to  experiments,  the  gas  manifold  and 
dosing  lines  were  heated,  under  vacuum,  to  100  °C  to  minimize 
water  contamination.  A  stainless  steel  reservoir  containing  30 
of  agent  was  then  connected  to  the  manifold,  and  the 
sample  was  purified  via  three  freeze— pump— thaw  cycles. 
Within  the  main  UHV  chamber,  controlled  exposure  of  the 
surface  to  the  agent  of  interest  was  achieved  through  a 
directional  doser  positioned  2  cm  from  the  center  of  the  silica 
sample.  Agent  uptake  was  monitored  via  IR  spectroscopy, 
performed  with  a  Bruker  Vertex  80v  FTIR  spectrometer 
coupled  to  the  UHV  chamber.  Each  spectrum  was  the  average 
of  100  scans  collected  at  a  2  cm“^  resolution  with  a  6  mm 
aperture  and  20  kHz  scanning  velocity.  The  spectra  presented 
here  are  difference  spectra,  where  the  reference  spectrum  is  that 
of  the  thermally  treated  silica. 

Previous  studies  have  shown  that  the  uptake  of  organo- 
phosphonate  simulants  on  the  surface  of  silica  occurs  through 
the  formation  of  hydrogen  bonds  to  free  silanol  groups.^’^ 
The  IR  spectral  signature  for  hydrogen  bonding  between 
adsorbates  and  silica  is  a  broad  SiO— H  stretching  mode  that 
appears  red-shifted  from  its  original  frequency.  Figure  2  shows 
IR  spectra  of  the  partially  dehydroxylated  Si02  particles  during 
adsorption  of  either  sarin  (Figure  2A)  or  soman  (Figure  2B). 
The  data  reveal  a  significant  change  in  the  original  free  SiO— H 
vibrational  band,  while  energies  of  other  modes  associated  with 
the  adsorbates  are  largely  independent  of  coverage.  These 
observations  indicate  that  adsorption  indeed  occurs  through 
strong  interactions  with  surface  hydroxyl  groups  and  that  the 
geometry  of  the  adsorbates  is  independent  of  coverage  within 
the  low-coverage  regime  employed  in  this  work. 
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Figure  2.  (A)  Difference  IR  spectra  of  amorphous  sihca  exposed  to 
gas-phase  sarin  (bottom)  and  IR  spectra  for  gas-phase  sarin  (top).  The 
spectra  displayed  in  the  lower  portion  of  the  graph  show  the  uptake  of 
sarin  during  exposure.  (B)  The  top  panel  shows  the  IR  spectram  for 
gas-phase  soman,  and  the  bottom  panel  shows  spectra  collected  during 
exposure  of  silica  to  soman.  The  IR  peaks  are  assigned  to  particular 
local  modes,  as  indicated  in  Table  1.  It  should  be  noted  that  definitive 
assignment  of  features,  especially  broad  peaks,  in  the  region  below 
~1400  cm~*  cannot  be  made  because  they  are  obscured  by  high 
absorbance  of  the  bulk  modes  of  silica  (see  Figure  6),  and  changes  in 
this  region  are  due  largely  to  the  change  in  index  of  the  refraction  of 
the  underlying  silica  upon  adsorption. 


The  IR  spectra  for  the  two  agents  are  very  similar,  reflecting 
the  similarities  in  their  molecular  structure  and  surface  bonding 
configurations.  Specifically,  Table  1  shows  that  the  low-energy 
bending  and  deformation  modes  are  within  ~10  cm”'  of  each 
other  for  the  two  adsorbates.  In  addition,  the  higher-energy  C— 
H  symmetric  and  asymmetric  modes  nearly  overlap  in 
wavenumber.  Finally,  the  primary  IR  modes  associated  with 
hydrogen-bonded  surface  hydroxyls  differ  by  only  ~40  cm”', 
which  suggests  that  sarin  and  soman  may  have  similar  binding 
structures,  energetics,  and  overall  mechanisms. 

For  both  adsorbates,  the  agent— surface  hydrogen  bond 
appears  to  be  dominated  by  interactions  between  the  surface 
SiOH  groups  and  the  sp^  oxygen  (P=0)  of  the  adsorbates. 
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Table  1.  Assignments  of  the  IR  Peaks  for  Sarin  and  Soman  in 
the  Gas  Phase  and  Adsorbed  on  Silica 


sarin  soman 


mode 

gas“ 

ads. 

gas“ 

ads. 

Si(0-H)fc. 

b 

3745 

b 

3745 

Si(0— H)bo„ded 

b 

3235 

b 

3196 

I/.fCHjC) 
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2994 
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t'.(CH3C) 

2949 
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2937 

p(CH) 
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2878 

2878 

Men,) 

b 

b 

1485 

1484 

Men,) 

b 

b 

1466 

1467 

Men,p) 
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1422 

1424 
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Men,c) 

1392 

1470 

b 

b 

Men,c) 

1387 

1461 

b 

b 

^(CHjC) 

1392 

1393 

b 

b 

4(CH3C) 

1387 

1381 

b 

b 

5.(CH3P) 

1326 

1325 

1326 

1324 

o 

II 

1303 

e 

1301 

e 

5.(CH3C) 

1387 

1461 

b 

b 

a(CH3C) 

1148 

1147 

i/(C-OP) 

1022 

1023 

1021 

1029 

CCC  bend 

991 

1000 

a(CH3) 

927 

931 

922 

927 

a(CH3) 

910 

910 

905 

905“^ 

Meec) 

887“^ 

886“^ 

869“* 

87U 

p(P-F) 

844 

846 

803 

825 

p(P-C) 

752 
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“See  ref  11.  ^’Indicates  that  modes  are  not  present  in  the  molecule. 
“^May  correspond  to  a  CCC  bend.  “^See  ref  12.  ‘^Obscured  by  bulk  Si02 
modes. 


This  conclusion  is  based  primarily  on  previous  theoretical 
studies,^  which  have  shown  that  the  largest  surface  binding 
energy  occurs  for  P=0"  H0Si  interactions.  For  example,  the 
hydrogen  bonding  energies  for  sarin  to  silanol  groups  were 
reported  to  be  approximately  4  kcal/mol  greater  when  bonding 
through  the  O(sp^)  than  those  through  the  O(sp^)  or  through 
the  F  groups  surrounding  the  molecule. 

Upon  hydrogen  bond  formation,  charge  transfer  from  the 
adsorbate  to  the  hydrogen  bond  acceptor  partially  populates 
the  a*  orbital  of  the  OH  group  to  weaken  this  bond.  This 
change  in  the  character  of  the  O— H  bond  is  reflected  in  the  red 
shift  of  the  vibrational  frequency  of  this  mode.  Furthermore,  a 
distribution  of  bonding  configurations  and  the  reduced  lifetime 
of  the  O— H  stretching  mode,  as  the  vibration  couples  to  other 
modes  of  the  adsorbate,  contribute  to  a  significant  broadening 
of  this  peak. 

Recent  work  has  suggested  that  the  shift  in  the  surface  SiO— 
H  stretching  frequency  is  strongly  correlated  to  the  strength  of 
the  hydrogen  bond  for  like  molecules  that  adsorb  to  silica 
through  a  central  phosphoryl  group.^  In  particular,  we  have 
shown  that  the  desorption  energies  from  silica  (for  the 
compounds  highlighted  in  Figure  l)  scale  linearly  with  the 
change  in  frequency  of  the  SiO— H  mode.^  Therefore,  one 
might  predict,  on  the  basis  of  the  IR  spectra  in  Figure  2,  that 
sarin  and  soman  have  similar  binding  energies.  We  have  tested 
this  hypothesis  through  quantitative  determination  of  the 
desorption  energy  for  these  two  molecules. 

The  experimental  temperature-programmed  desorption 
(TPD)  data,  which  display  the  rate  of  desorption  for  sarin 
and  soman  as  a  function  of  surface  temperature,  are  shown  in 
Figure  3A  and  C.  In  these  studies,  a  doubly  differentially 


Temperature  (K)  Coverage  (ML) 


Figure  3.  (A)  TPD  data  collected  for  the  desorption  of  sarin  from  amorphous  silica.  After  ejqjosure  of  the  silica  to  gas-phase  sarin,  the  surface 
temperature  was  ramped  at  a  rate  of  0.5  K/s.  (B)  The  desorption  energy  of  sarin  from  silica  as  a  function  of  surface  coverage  determined  by  inversion 
of  the  Polanyi— Wigner  equation.  (C)  TPD  data  for  soman  desorbing  from  amorphous  silica.  The  data  in  panel  C  was  analyzed  via  inversion  of  the 
Polanyi— Wigner  equation  to  provide  the  coverage-dependent  desorption  energy  (D). 
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pumped  mass  spectrometer  was  used  to  monitor  the  desorption 
rate  for  each  molecule  during  a  linear  thermal  ramp  of  the 
surface.  We  note  that  the  mass  spectrometric  cracking  patterns 
for  sarin  and  soman  molecules  as  they  desorbed  from  the 
surface  were  identical  to  those  of  the  parent  molecules,  which 
indicates  that  molecular  dissociation  is  very  limited  on  the  silica 
surface.  This  conclusion  is  supported  by  the  IR  data  of  Figure  2, 
which  shows  spectroscopic  evidence  for  molecular  adsorption. 

Determination  of  the  desorption  energy  (£j)  from  the  TPD 
data  was  achieved  through  inversion  of  the  Polanyi— Wigner 
equation 


dt 


5,  rj 


(1) 


In  this  equation,  u  is  the  preexponential  factor,  ©  is  the 
coverage,  is  the  surface  temperature,  fcg  is  the  Boltzmann 
constant,  and  n  is  the  order  for  desorption  (simple  molecular 
desorption  is  expected  to  be  first  order,  n  =  l).^^  The  inversion 
analysis  is  accomplished  by  solving  the  above  equation  for  the 
coverage-dependent  desorption  energy 


Eii@)  =  -hT^  In 


(- 


d©/dt\ 
iy@  j 


(2) 


For  this  work,  we  used  a  pre-exponential  factor  determined  for 
the  desorption  of  the  simulant  DMMP  from  the  silica  sample, 

=  4.0  X  s~h^  Although  safety  constraints  precluded  us 

from  performing  the  types  of  repeated  exposure-dependent 
TPD  studies  required  for  determination  of  a  unique  pre¬ 
exponential  factor  for  the  agents,  the  pre-exponential  for 
DMMP  has  been  shown  to  accurately  model  the  TPD  data  for 
all  of  the  simulants  highlighted  in  Figure  l7  Close  agreement  of 
the  agent  data  with  previously  observed  trends  in  the 
desorption  energy  for  the  simulants,  as  described  below, 
helps  to  justify  our  use  of  the  pre-exponential  factor  for  DMMP 
in  the  current  TPD  analysis. 

The  Ej(ff)  results  (Figure  3B  and  D)  from  inversion  of  the 
TPD  data  produce  a  range  of  binding  energies,  well  above  45 
kj/mol,  for  both  agents.  In  fact,  the  desorption  energy  for 
soman  (blue)  approaches  60  kJ/mol  at  the  lowest  coverages. 
The  high  activation  energies  of  desorption  found  in  this  work 
exceed  typical  hydrogen  bond  strengths;  however,  the 
energetics  are  consistent  with  recent  work  by  a  number  of 
groups  that  have  performed  electronic  structure  calculations  of 
many  agents  and  simulants  on  the  surface  of  silica.*'^'^^  They 
are  also  consistent  with  calorimetric  measurements  of  the  heats 
of  adsorption  for  sarin  on  silica.*  The  high  hydrogen  bond 
energy  appears  to  be  driven  by  the  significant  electron  density 
surrounding  the  sp^  oxygen  atom  of  the  agents.  Even  with  the 
strongly  electron  withdrawing  fluorine  substituent  on  the 
molecules,  the  charge  on  the  oxygen  may  be  as  high  as  —0.6  e, 
where  e  is  the  charge  of  1  electron,  as  determined  from 
electronic  structure  calculations.^  The  significant  charge  on  the 
oxygen  creates  large  electrostatic  interactions  with  the  polar 
surface  hydroxyl  groups  and  facilitates  partial  charge  transfer  (as 
evidenced  by  the  large  red  shift  of  the  SiO— H  stretch)  to 
produce  strong  hydrogen  bonding  forces  (as  evidenced  by  the 
TPD  measurements). 

Overall,  the  data  of  Figure  3  show  that  the  molecule— surface 
interactions  are  weakest  for  the  highest  surface  coverages  (~45 
kJ/mol  for  sarin,  ~47  kJ/mol  for  soman)  but  become  stronger, 
by  ~5  kJ/mol,  as  the  coverage  decreases.  This  observation 
indicates  that,  under  these  experimental  conditions,  the 


molecules  readily  diffuse  to  the  most  stable  unoccupied  sites 
that  produce  the  tightest  binding.  In  this  way,  the  increase  in 
energy  suggests  that  there  are  a  variety  of  surface  sites  or 
configurations,  each  leading  to  slightly  different  energies. 
Alternatively,  there  may  be  small  repulsive  interactions  between 
the  adsorbates,  the  extent  of  which  depends  on  coverage.  The 
strongest  binding  sites,  responsible  for  the  upturn  in  the  Ej 
curve  for  soman  at  low  coverage,  may  be  defect  sites  or  areas 
where  many  hydroxyl  groups  contribute  simultaneously  to  the 
adsorption.  The  binding  energy  for  isolated  hydroxyl  groups  in 
defect-free  regions  of  the  surface  can  be  estimated  by 
extrapolating  the  linear  portion  of  the  energy  distribution 
curves  to  zero  coverage  (dashed  line  in  Figures  3B  and  D).  This 
analysis  indicates  that  the  energy  for  desorption  of  sarin  from 
isolated  SiOH  groups  is  50  ±  2  kJ/mol  and  that  of  soman  is  52 
±  2  kj/mol. 

Recently,  Wilmsmeyer  et  al.^'^  explored  the  relationship 
between  molecular  functionality  and  interfacial  organophosph- 
onate  hydrogen  bonding  at  the  gas— surface  interface  of  sihca. 
In  their  studies,  the  coverage-dependent  desorption  energies  for 
the  CWA  simulants  DMMP,  DIMP,  DMCP,  MDCP,  and  TMP 
were  measured  in  conjunction  with  IR  spectroscopic  studies  of 
the  adsorbed  species.^’^  The  reported  Ej  values  were  43.0  ±  0.8 
(MDCP),  48.4  ±  1.0  (DMCP),  52.4  ±  0.6  (TMP),  54.5  ±  0.5 
(DMMP),  and  57.9  ±  1  (DIMP)  kj/mol.  The  experimental 
measurements,  together  with  electronic  structure  calculations, 
revealed  that  the  desorption  energy  from  an  isolated  silanol 
group  is  directly  related  to  the  charge  on  the  sp^  oxygen  atom 
of  the  simulant,  thus  demonstrating  that  hydrogen  bond 
formation  occurred  exclusively  through  the  central  phosphoryl 
group.  The  other  substituents  affect  the  charge  on  the  sp^ 
oxygen,  but  they  are  not  directly  involved  in  hydrogen  bond 
formation.  That  is,  the  binding  energy  decreases  with  the 
electron-withdrawing  ability  of  the  surrounding  functional 
groups.  For  example,  the  chlorine  substituents  of  MDCP 
significantly  decrease  the  electron  density  of  the  phosphoryl 
moiety,  thus  weakening  the  interfacial  hydrogen  bond. 
Substitution  of  the  chlorine  with  a  methyl  or  large  isopropoxy 
group  (e.g.,  DIMP)  results  in  stronger  binding  with  the  surface. 
Similar  results  were  reported  in  the  recent  computational 
studies  by  Troya  et  al.,^  in  which  they  demonstrated  that  the 
fluorine  substituent  in  sarin  reduced  the  electron  density  of  the 
sp^  oxygen  but  that  this  effect  was  largely  countered  by  the 
isopropoxy  group  of  the  molecule,  which  accounted  for 
computed  binding  energies  similar  to  those  of  the  most 
common  simulants.  The  net  effect  of  the  halogen  atoms 
appears  to  be  that  of  withdrawing  electron  density  from  the  sp^ 
oxygen,  which  is  the  most  important  factor  in  determining  the 
strength  of  adsorption  and  the  activation  energy  for  desorption. 

The  relationship  between  the  energy  for  desorption  and  the 
change  in  the  vibrational  frequency  of  the  surface  SiO— H  mode 
upon  adsorption  has  been  shown  to  be  remarkably  linear  for 
the  series  of  simulants  presented  in  Figure  1.  With  the  results 
reported  above,  we  find  that  the  actual  agents  agree  very  well 
with  this  trend.  Figure  4  shows  the  previously  published 
relationship  between  desorption  energy  and  SiO— H  vibrational 
frequency  for  the  simulants,^’ along  with  the  addition  of  two 
key  data  points  for  the  agents.  This  relationship  demonstrates 
that  the  binding  mechanism  for  each  molecule  is  similar,  with 
the  dominant  contribution  to  the  strength  of  the  hydrogen 
bond  being  interactions  between  the  sp^  oxygen  of  the 
adsorbates  and  the  SiOH  groups  of  the  surface.  Further,  the 
simulant  that  most  closely  mimics  the  energetics  of  the  agents 
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Figure  4.  Plot  of  the  change  in  wavenumber  of  the  SiO— H  stretching 
mode  (At/)  and  the  desorption  energy  (£j)  for  the  agents  and 
simulants. 


is  TMP.  Interestingly,  the  strong  electron-withdrawing 
character  of  the  fluorine  atom  on  the  agents  does  not  lead  to 
weak  hydrogen-bonding  interactions  because  this  effect  appears 
to  be  balanced  by  the  charge-donating  character  of  the  oxygen- 
containing  substituents. 

This  correlation  between  nerve  agents  and  the  most  common 
simulants  will  better  enable  researchers  to  choose  the  most 
appropriate  simulants  as  well  as  provide  critical  benchmarks  for 
computational  studies  of  agents  where  hydrogen-bonding 
interactions  play  a  role  in  the  overall  chemistry.  Most 
importantly,  the  results  provide  fundamental  insight  into  the 
structure— function  relationship  in  hydrogen-bond  formation  at 
the  gas— surface  interface  that  extends  our  understanding  of  this 
critical  chemical  interaction. 

■  EXPERIMENTAL  METHODS 

The  custom  multifunctional  UHV  chamber  specifically 
designed  to  safely  allow  experiments  with  ultratoxic  chemicals 
used  for  this  work  is  described  in  a  previous  article^  and  is 
depicted  schematically  in  Figure  5.  The  chamber  combines 
multiple  vacuum  dosing  methods  for  liquids  and  vapors  and  the 
ability  to  perform  surface-sensitive  IR  spectroscopy  and  XPS  as 
well  as  mass  spectrometry  (MS).  The  instrument  is  built 
around  a  23  L  UHV  chamber  constructed  out  of  316L  stainless 
steel  with  all  ports  equipped  with  con-flat  flanges  (Kurt  J. 
Lesker  Company).  The  chamber  is  evacuated  with  a  2200  L/s 
magnetically  levitated  turbomolecular  pump  (A2200C,  Edwards 
Vacuum).  This  pump  is  backed  by  a  scroll  pump  (XDS-35, 
Edwards  Vacuum). 

Samples  are  introduced  into  the  chamber  via  a  high-vacuum 
load-lock  system  that  is  coupled  to  the  chamber  through  a  6  in. 
gate  valve.  Samples  are  mounted  onto  a  molybdenum  platen, 
which  is  in  turn  secured  onto  the  linear  transfer  arm  via  a 
spring-loaded  fork  (customized  STLC  transfer  system, 
Thermionics,  Inc.).  Once  evacuated,  the  load-lock  gate  valve 
is  opened  to  enable  transfer  of  the  sample  platen  on  to  the 
manipulator  (Thermionics,  Inc.),  which  allows  travel  in  the  X, 
Y,  and  Z  dimensions,  full  360°  of  rotation,  and  tilt  in  the  X  and 
Y  planes.  This  enables  the  sample  to  be  positioned  as  required 
for  the  experiment.  The  sample  platen  features  a  transferrable 
type  K  thermocouple,  which  is  spot-welded  to  the  surface. 

Sarin  and  soman  were  dosed  as  a  neat  vapor  using  a 
multivalved  stainless  steel  high-vacuum  manifold.  The  entire 
path  of  the  manifold  was  evacuated  and  heated  to  100  °C  prior 
to  each  experiment  to  eliminate  water  and  other  contaminants. 


Figure  5.  Schematic  (not  to  scale)  of  the  UHV  system  showing  the 
primary  components:  (A)  sample  transfer  translation  stage,  (B)  load- 
lock  chamber  and  pump,  (C)  chemical  agent  vapor  dose  manifold,  (D) 
gate  valve,  (E)  IR  hght  entrance  optics,  (F)  capillary  array  doser,  (G) 
apertures  and  differential  pumping  stage,  (H)  mass  spectrometer 
chamber  and  pump,  (l)  cryostat  for  sohd  sorbent  doser,  (j)  directional 
agent  vapor  doser,  (K)  IR  light  exit  optics  and  detector,  (L)  XPS 
system  and  ion  gun,  (M)  forehne  pumps  for  the  turbomolecular 
pumps. 


The  vapor  pressure  of  the  CWAs  was  sufficient  to  allow  room- 
temperature  vapor  dosing  once  the  entire  manifold  volume  was 
evacuated.  Freeze- pump— thaw  cycles  using  liquid  nitrogen 
removed  dissolved  gases  as  well  as  remaining  air  from  the 
manifold  prior  to  vapor  dosing.  The  vapor  manifold  line 
introduced  CWA  vapor  to  the  UHV  chamber  through  a 
precision  leak  valve  with  a  directional  heated  tube  on  the  UHV 
side  to  direct  the  vapor  to  the  sample  surface. 

Uptake  of  the  CWA  was  monitored  in  situ  using  RAIRS.  Eor 
these  data,  a  vacuum  IR  spectrometer  (Vertex  80v,  Bruker)  has 
been  modified  to  mate  with  the  UHV  system.  Specifically, 
radiation  (SiC  glowbar)  exits  the  spectrometer  and  is  directed 
to  a  custom  external  optics  compartment  that  includes  a/=  250 
mm  parabolic  mirror  for  directing  the  beam  through  a 
differentially  pumped  KBr  window  and  focusing  it  onto  the 
surface  sample.  The  beam  size  at  the  center  of  the  chamber  is 
estimated  to  be  approximately  9  mm  X  6  mm.  The  angle  of 
incidence  for  reflection  from  a  flat  surface  is  86°  with  respect  to 
the  sample  plane.  After  exiting  the  chamber  through  a  second 
differentially  pumped  KBr  window,  the  diverging  reflected 
beam  is  then  focused  with  an  ellipsoidal  mirror  (/1//2  =  250/40 
mm)  onto  the  sensor  element  of  a  liquid-nitrogen-cooled 
mercury  cadmium  telluride  (MCT)  detector.  The  entire  beam 
path  can  be  operated  under  vacuum,  which  is  important  for 
removing  background  gases  from  the  spectra.  In  addition,  the 
mirrors  in  the  system  are  high  reflectivity  Au-coated  mirrors  to 
minimize  loss  of  photons  throughout  the  beam  path. 

The  IR  spectrum  of  silica  is  shown  in  Figure  6,  which  was 
collected  after  pretreatment  of  the  sihca  at  700  K  for  5  min. 
The  spectrum  shows  a  sharp  peak  at  3745  cm”*,  which  is 
characteristic  of  isolated  hydroxyls.  The  absence  of  a  broad 
feature  near  3500  cm”*  indicates  the  removal  of  most 
associated  hydroxyls.  Work  by  Zhuravlev^  suggests  that, 
under  our  annealing  conditions,  free  hydroxyls  dominate  the 
surface,  but  there  remain  approximately  20%  geminal  groups. 
The  small  peak  near  3000  cm”*  wavenumbers  is  likely  due  to 
trace  hydrocarbon  contamination  that  could  not  be  removed 
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Figure  6.  IR  spectrum  of  silica  after  pretreatment  at  700  K  for  5  min. 


regardless  of  sample  annealing  conditions.  The  intensity  and 
position  of  this  peak  were  unaffected  by  sarin  and  soman 
adsorption. 

Although  IR  spectroscopy  yields  information  about  the 
vibrational  frequencies  of  surface  adsorbates,  complete 
characterization  of  the  activation  energy  for  desorption  was 
performed  via  kinetic  analysis  of  TPD  data.  TPD  is  performed 
by  first  exposing  the  surface  to  a  species  of  interest.  During  this 
initial  exposure,  RAIRS  is  employed  to  monitor  surface  uptake 
and  determine  coverage.  Initial  dosing  is  typically  performed  at 
low  surface  temperatures  where  residence  times  are  sufficient  to 
ensure  constant  coverage  following  dosing.  Once  dosing  is 
complete,  the  chamber  is  evacuated  to  remove  background  gas- 
phase  agent.  The  temperature  of  the  surface  is  then  increased 
by  driving  the  sample  heater  with  a  power  supply  under 
proportional-integral-derivative  control  (LPS-800-1,  Ther- 
mionics).  With  the  liquid  nitrogen  reservoir  on  the  sample 
holder  filled,  the  temperature  can  be  ramped  at  a  rate  of  0.5  K/s 
with  approximately  ±1  K  precision.  During  the  thermal  ramp, 
IR  and  mass  spectra  are  recorded  to  track  desorption  from  the 
surface.  Mass  spectral  intensity  versus  surface  temperature 
represents  the  rate  of  desorption,  which  is  governed  by  the 
activation  energy  for  desorption.  In  the  absence  of  chemical 
reactions,  the  desorption  rate  curves  can  be  analyzed  to  reveal 
properties  about  the  activation  energy  for  desorption. 

As  products  are  released  from  the  surface  during  TPD,  they 
are  tracked  with  a  doubly  differentially  pumped  quadrupole 
mass  spectrometer,  (Extrel,  MAXIOOO  m/z  =  1— 1000).  Each 
differential  stage  is  evacuated  by  a  450  L/s  magnetically 
levitated  turbomolecular  pump  (Edwards  STP-451C),  which 
enables  up  to  3  orders  of  magnitude  pressure  differential 
between  the  main  chamber  and  mass  spectrometer  chamber 
during  an  experiment.  The  line-of-sight  of  the  mass 
spectrometer  is  defined  by  a  series  of  apertures  that  separate 
the  differential  pumping  stages  such  that  the  ionizer  of  the  mass 
spectrometer  views  a  1  cm^  spot  on  the  surface  when  the 
surface  is  located  at  the  focal  point  of  the  main  chamber.  The 
acceptance  angle  of  the  spectrometer  is  1°  in  polar  and 
azimuthal  angles.  As  indicated  in  Figure  5,  the  mass 
spectrometer  is  aligned  such  that  it  views  the  surface  during 
dosing  and  while  performing  IR  spectroscopic  measurements. 
In  this  way,  changes  in  surface  vibrational  modes  and 
concentrations  during  dosing  can  be  directly  correlated  with 
the  mass  spectra  of  species  desorbing  from  the  surface  at  any 
point  in  time.  A  second  mass  spectrometer  (RGA  300M, 
Stanford  Research  Systems,  m/z  =  1—300)  mounted  directly  on 


the  UHV  chamber  provided  independent  confirmation  of  intact 
molecular  desorption  of  the  CWA  and  was  used  to  confirm 
high  purity  of  the  agent  as  the  vapor  entered  the  chamber. 

The  absence  of  significant  CWA  decomposition  on  the  silica 
surfaces  during  the  experiments  was  corroborated  by  XPS  data 
recorded  following  the  IR/TPD  experiments.  The  XPS 
instrumental  component  of  this  UHV  system  consists  of  a 
dual  anode  (A1  Ka  and  Mg  Ka)  X-ray  source  (DAR  400, 
Omicron  Nanotechnology)  and  is  equipped  with  a  Z 
translation  stage  with  5°  X—Y  tUt  for  alignment  with  the 
sample.  Detection  of  photoelectrons  is  achieved  with  an 
Omicron  Sphera  hemispherical  analyzer  with  a  five-channel 
detector  mounted  at  a  45°  angle  relative  to  the  X-ray  source. 
Energy  calibration  was  performed  using  the  Au  4f7/2  peak 
observed  from  the  gold  substrate  at  84.00  eV.  A  typical  XPS 
data  set  is  shown  in  Figure  7,  which  was  collected  from  a  silica 


Figure  7.  Survey  XPS  data  collected  after  several  cycles  of  sarin  and 
soman  exposure  followed  by  TPD.  Note  the  absence  of  P  or  F  signal. 
Au  is  present  due  to  the  substrate  supporting  the  sdica  film,  and  the 
small  amount  of  carbon  is  due  to  trace  concentrations  of  hydrocarbons 
on  the  sihca  surface. 


surface  after  several  TPD  cycles  using  both  sarin  and  soman. 
The  absence  of  phosphorus  or  fluorine,  which  would  be 
expected  to  be  present  on  the  surface  of  metal  oxides  after 
reacting  with  organophosphorus  compounds,^^*^’^*  confirms  that 
desorption  is  molecular  and  does  not  involve  significant 
chemical  reactions. 
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